1 



(N 
O 
O 
(N 

< 

(N 



> 
O 

O 

oo 

o 

O 



X 



Large N reduction with overlap fermions 

J. Kiskis a , R. Narayanan b , and H. Neuberger *''' 
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b Department of Physics, Florida International University, Miami, FL 33199 
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We revisit quenched reduction with fermions and explain how some old problems can be avoided using the 
overlap Dirac operator. 



1. Introduction 

This contribution, based on joint work, pro- 
vides some background to Narayanan's talk 

Lattice field theory provides an approach to 
calculating numerical values or bounds for physi- 
cal quantities from first principles without relying 
on perturbation theory or on non-systematic ap- 
proximations. However, most of the numerical 
results about strong interactions involve fermions 
and are obtained in valence (quenched) QCD. 
This is not a systematic approximation because 
the zero flavor limit is singular. As a result, many 
lattice results still do not have a status higher 
than estimates coming from other phenomenolog- 
ical approximations that make ad-hoc approxima- 
tions of different types. It makes little sense to go 
beyond few percent accuracy in simulations with 
valence QCD. 

The long term goal of our work is to replace 
valence QCD with planar QCD in either the 't 
Hooft (|] or the Veneziano || limit. Unlike va- 
lence QCD, for many quantities of interest, pla- 
nar QCD is believed to be the first term in a sys- 
temtic series. Our initial focus is on the two point 
functions of quark bilinears. In the planar limit, 
such a correlation function is characterized by an 
infinite number of pole masses and correspond- 
ing residues. Using the lattice, it should be pos- 
sible to obtain good values for these correlation 
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functions in momentum space at Euclidean mo- 
menta that are smaller than the inverse lattice 
spacing. Such information would be phenomeno- 
logically useful, for example, in the analysis of 
weak semi- and non-leptonic decays pj] . 

2. The planar limit 

Let us consider QCD with gauge group SU (N c ) 
instead of SU(3). We are in Euclidean space, 
and assume that the volume V is already taken 
to infinity. Similarly, the ultraviolet cutoff A has 
also been eliminated. The coupling constant has 
disappeared too, generating an intrinsic scale in 
terms of which all predictions are pure numbers 
depending only on N c . There are good choices 
for this intrinsic scale, in the sense that they lead 
to simple N c dependences. Feynman diagrams 
and more general counting arguments tell us what 
powers of N c to scale out so that remaining quan- 
tities have finite limits at N c = oo. In short, a 
nonperturbatively defined large N c limit seems to 
exist and is reachable by lattice techniques. 

In perturbation theory, the large N c limit is 
easiest to understand before renormalization is 
taken into account. To get planar perturbative 
QCD one scales the gauge coupling constant g 2 
by a factor of N c and N c is taken to infinity keep- 
ing A = g 2 N c fixed. Using the known terms in 
the /3-function one easily derives how a perturba- 
tive scale can be defined so that it stays finite at 
infinite N c . 

The above 't Hooft limit generalizes to the case 
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where the number of flavors, Nf also is taken to 
infinity with the ratio P = jf kept fixed. This 
gives the Veneziano limit. 

On the lattice, fast convergence to the N c = oo 
limits was observed for the string tension, glucball 
mass, topological susceptibility, and finite tem- 
perature transition point in pure gauge theory . 
These calculations are of theoretical interest; for 
practical applications it is certainly cheaper to 
just work at N c = 3. 

Having learned that the large N c limit is mean- 
ingful outside perturbation theory and that it 
provides numbers quite close to N c — 3, we note 
that large N c reduction Q provides a potential 
shortcut that would enable us to get the pla- 
nar limit numbers in a substantially more effi- 
cient way because of the dramatic reduction in 
the number of degrees of freedom one needs to 
integrate over. In particular, the fermionic ma- 
trix becomes so small that it can be stored in its 
entirety in memory and dealt with in a more di- 
rect way than in standard simulations. 

3. Quenched Reduction 

Reduction works on lattices with finite num- 
bers of sites. In physical applications, reduction 
implies a change in the order of limits of the tradi- 
tional jj- expansion. Now, we first take N c (pos- 
sibly also Nf) to infinity and only subsequently 
deal with the large V and large A limits. 

The basic content of reduction is that there 
exists a certain "prototype" lattice gauge the- 
ory which encompasses lattice gauge theories with 
varying group sizes on differently sized lattices. 
However, the entire prototype family has the 
same planar limit. Thus, large V can be traded 
for larger N c . The main point is that this tradeoff 
appears to be very lucrative. 

The prototype variables are d S'C7(7V)-matrices 
T M , /u = 1, 2.., d and two fermionic N x M com- 
plex matrices "J and "P. The pure gauge action is 
given by 

Sg = Tr {C^ U C^ V ) , C M „ = [T^, T v \. (1) 

It is invariant under SU(N) conjugation T M — » 
fiT Al fi t and under a Z(N) d acting by -> 



To add fermions one introduces a "Dirac" oper- 
ator D(T), which also transforms by conjugation, 
and a term S f = Tr (*D(T)*) This breaks the 
Z(N) d . There is an additional fermionic symme- 
try under the obvious action of SU(M). 

Assuming that N and M factorize, N — 
nLiL2---L c i and M = mL\L2-..L^ one can con- 
strain the basic variables in the prototype so that 
one obtains a traditional lattice gauge theory on 
a toroidal Lj x L 2 - x lattice. A planar limit 
can be defined on the lattice by scaling the lat- 
tice gauge coupling with n in the 't Hooft fashion. 
The basic claim of [|| was that for a certain class 
of interesting observables, the differently con- 
stained prototypes become identical in this pla- 
nar limit, at finite lattice spacing. However when 
the continuum limit is approached, the identity 
is spoiled by a phase transition. Quenched re- 
duction Jjj] is designed to eliminate this effect by 
applying another constraint to the prototype the- 
ory: The eigenvalues of the T„ matrices are frozen 
at d independently drawn sets of uniformly dis- 
tributed points on the unit circle parametrized by 
angles Q^. Further ||, the T M matrices entering 
D(T) are multiplied by U(l) phase factors of the 
form e lPfl . Observables obtained by carrying out 
a quenched average with fixed angles have to be 
subsequently averaged over uniformly distributed 
9^ and p^. A large N limit can be defined as 
before, scaling the lattice coupling with N. 

4. The resolution of some old problems 

In the past, when chirality and the lattice were 
thought to be irreconcilable, the preferred way 
to introduce fermions was in a continuum version 
of quenched reduction ||] . In this version one re- 
places the unitary matrices of the prototype by 
hcrmitian matrices: 7), — » 1 + a^, the eigen- 
values of , are frozen and drawn from a uniform 
distribution over an interval [— py]. Similarly, 
the Dirac operator D(T) gets replaced by a D(A) 
and the random variable p^ is introduced by re- 
placing A by A^ + with randomly drawn 
again from [— A, A]. 

In perturbation theory this looks fine. Surpris- 
ingly, there seems to be no clash between the UV 
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regularization and gauge invariance even when 
the continuum theory would be chiral and anoma- 
lous. In addition, there seems to be no room in 
even d = 2k for topological charge, since the nat- 
ural object 



^fJ,lVlfJ,2V2----fJ.kVk 



Tr ( , A Vl ] [A^ 2 , A U2 ].... [A flk , A Vk ] ) 



(2) 



vanishes identically. 

In the past, it was unclear how this regulated 
continuum reduced model differed from the lat- 
tice version. It made some sense that topology 
would become a strange concept because the dis- 
appearance of spacetime took away the concept 
of a local topological density. Also, so long as 
fermions were dealt with in the 't Hooft limit, the 
issue of anomalies did not look pressing. Never- 
theless, there was enough uncertainty surround- 
ing this issue that no further work on fermions in 
reduced models was done. The lattice provided 
no attractive alternative. 

The main new observation in this context is 
that choosing D(T) as the overlap Dirac oper- 
ator [jl0| frees the lattice version with fermions 
from all difficulties. The overlap also provides a 
definition of topological charge in the pure 
gauge reduced model, opening the possibility to 
be consistent with recent studies of topology on 
the lattice at large N c in traditional formula- 
tions. Moreover, topological obstructions that 
reflect anomalies are known to exist in ordinary 
lattice gauge theories already on finite lattices, as 
emphasized in |12] ]. These obstructions will also 
appear in the reduced model. Thus, so long as 
we stay on the lattice, quenched reduction will 
not eradicate anomalies. 

5. Comments 

Reduction seems to be closely related to the 
technique of "field orbifolding" ]l3|] . For a pure 
U (N c ) gauge theory, the constraints that turn the 
prototype model into a regular lattice gauge the- 
ory are equivalent to imposing invariance under 
a subgroup of the symmetry group of the proto- 
type. It is known that field orbifolding preserves 
the planar limit in certain continuum filed theo- 
ries and has a geometrical interpretation in string 



theory. Could one gain new insights into reduc- 
tion from this? 

Numerically, twisted reduction jl4| used to 
be considered superior to quenching in the pure 
gauge sector. Unfortunately, the introduction of 
fermions is problematic |^5| . Can one get around 
the fermion problem in some way while sticking to 
twisted reduction instead of quenched reduction 
in the pure gauge sector? 

6. Conclusions 

This is a good time to numerically revisit 
quenched reduction. Some of the old problems 
have been solved: Instanton effects and anomalies 
can be incorporated. Presently available compu- 
tational power, albeit still short of what is needed 
for a full QCD simulation, might be amply capa- 
ble of dealing with reduced models. 
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